arXiv:1507.05626v2 [astro-ph.GA] 22Jul2015 


Mon. Not. R. Astron. Soc. 000, [^?? (2002) 


Printed 23 July 2015 (MN style file v2.2) 


THE MAGELLANIC INTER-CLOUD PROJECT 
(MAGIC) II: SLICING UP THE BRIDGE 

N. E. D. NoeO*, B. C. Conn^, J. I. Read^, R. Carrera^’^, A. Dolphin^, and H.-W. Rix® 

^Department of Physics, University of Surrey, Guildford, GU2 7XH, UK 
^Gemini Observatory, Recinto AURA, Golina El Pino s/n. La Serena, Ghile 
^Institute de Astrofisica de Ganarias, C/ Via Ldctea s/n, 38200, La Laguna, Tenerife, Spain 
Departamento de Astrofisica, Universidad de La Laguna, Tenerife, Spain 
^Raytheon Company, PO Box 11337, Tucson, AZ, 85734-1337, United States of America 
^ Max Planck Institut fiir Astronomic, Konigstuhl 17, 69117, Heidelberg, Germany 


Released 2002 Xxxxx XX 


ABSTRACT 

The origin of the gas in between the Magellanic Clouds (MCs), known as the 
Magellanic Bridge (MB), has always been the subject of controversy. To shed light 
into this, we present the results from the MAGIC II project aimed at probing the 
stellar populations in ten large fields located perpendicular to the main ridge-line of 
HI in the Inter-Cloud region. We secured these observations of the stellar populations 
in between the MCs using the WFI camera on the 2.2 m telescope in La Silla. Using 
colour-magnitude diagrams (CMDs), we trace stellar populations across the Inter- 
Cloud region. In good agreement with MAGIC I, we find significant intermediate- 
age stars in the Inter-Cloud region as well as young stars of a similar age to the 
last pericentre passage in between the MCs (^200 Myr ago). We show here that the 
young, intermediate-age and old stars have distinct spatial distributions. The young 
stars correlate well with the HI gas suggesting that they were either recently stripped 
from the SMC or formed in-situ. The bulk of intermediate-age stars are located mainly 
in the bridge region where the HI column density is higher, but they are more spread 
out than the young stars. They have very similar properties to stars located ~ 2 Kpc 
from the SMC centre, suggesting that they were tidally stripped from this region. 
Finally, the old stars extend to some 8 Kpc from the SMC supporting the idea that 
all galaxies have a large extended metal poor stellar halo. 

Key words: Galaxies: evolution, galaxies: stellar populations, stars: evolution. 


1 INTRODUCTION 

A key challenge in modern astrophysics is to understand 
how galaxies form and evolve. In particular, interactions be¬ 
tween galaxies are an important driver of galaxy formation 
and evolution (e.g. Larson fc Tinsley|1978 Barnes & Hern 
|quist|[T992[ ). These encounters reshape galaxies by transfer¬ 
ring mass, energy, gas and metals between them, creating 
new sites of star formation both within the galaxies and in 
stripped gas amidst them. Important clues about the galaxy 
formation process lie in the faint edges of galaxies. Galaxy 
formation simulations suggest that almost all galaxies, even 
the smallest, should contain an extended old, metal poor, 
stellar halo that is primarily composed of the detritus from 


galactic mergers (e.g. White & Rees 1978 Searle & Zinn 
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1978[ [Bullock fc Johnston|2005|[Read et al.|2006|[Uooper et| 


al.pOlO I. This is supported by a host of observations from 

large spirals (e.g. Roman 1954| Eggen et al. 1962 [Kalirai 


e£^JJ200^ McConnachie A al.pOOOTChnpmA^et^^ajj2006| 


Gilbert et al. 2012[ Trujillo fc Bakos| 20131 Monachesi et| 


al. 2013| Deason et al. 2014| Greggio et al. [2014 1 to dwarf 


galaxies (e.g. Aparicio et al.pOOO Minniti et al.|20^ Noel 


fc Gallart|2067 Pejcha fc Stanek|2009 Strader et al. [2012 1. 
Such halos retain a fossil record of past mergers and inter¬ 


actions in the form of extended tidal debris (e.g. Searle fc 


^ Martmez-Delgado et al.|2001[ [Munoz et al.|200^ 

Sohn et al.j2007} ' Amorisco et al.||2014 K 


The distribution and extent of such debris, or a char¬ 
acteristic underlying old stellar halo are both sensitive to 


the total mass and size of a galaxy (e.g. Noel fc Gallart 


|2007| [Noel et al.||2013a] [Deason et al.|[^14| ) . Thus, probing 

the faint edges of galaxies can also give us constraints on 
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their dark matter content. Observing these remote edges of 
galaxies is a complicated task due to the faintness of these 
outlying regions, typically /i^28 mag/arcsec^ (e.g. Noel & 


Gallart|[2007 1. Reliable surface brightness measurements at 


these faint levels require accurate flat held division, sky de¬ 


termination, bright star masking, among others (e.g. Barker 
et al.|2009| ). 

Our nearest irregular neighbours, the Magellanic 
Clouds (MCs) constitute a unique opportunity to study 
a gas rich interacting system close-up and to investigate 
faint outskirts of galaxies. First, their close proximity (~50- 
70Kpc away) enables us to resolve their individual stars, 
which allows to directly probe fainter surface brightness lev¬ 
els placing constraints on ages and metallicities without the 


need to use stellar population synthesis (e.g. Noel et al. 


2013b I. Second, due to their mutual association over the 


past few Gyrs (as determined from their proper motions: 


Kallivayalil et al. 

20061 Besla et al. 

2007 

IPiatek et al. 

2008 

^ost^et""aLll20TT 

IBesla et al.|2012 

iKal 

ivayalil et al. 

2013 

van der Marel & 

Callivayalil|2014[ Cioni et al.|20141 and nu- 

merical simulations: Bekki [20111 [Diaz & Bekki||2012| [Besla 


(e.g. Stanimirovic et al.|2000| Kemper et al.|2010|| Bolatto et 


al.|2011 Matsuura et al.||2013[ ). Indeed various gaseous fea¬ 

tures were created by the interactions of the MCs with each 
other and with the Milky Way (MW) such as the Magellanic 


Bridge (MB; Hindman et al. 19631, a common envelope of 
HI spanning ~13 Kpc in which both the Large Magellanic 
Cloud (LMC) and the Small Magellanic Cloud (SMC) are 
embedded; the Magellanic Stream, a pure gaseous stream 


trailing the MCs as they orbit the MW ( Mathewson et al. 
1974 Nidever et al. 12010 1 connecting both of the MCs with 
the MW, and the Leading Arm, a mass of high-velocity 
clouds leading the MCs | Putman||2000 K 

Since its discovery, the MB was believed to have a 
tidal origin, forming during an encounter between the MCs 
I Mathewson|1985 1. The SMC Wing -the Eastern protuber¬ 
ance of the SMC main body- extends towards the LMC into 
the Bridge ( Irwin et al.||1990 ) indicating that this tidal in¬ 
teraction could have formed a disk around the SMC that 
was later on stripped to form the MB | Grondin et al.|1992 1. 
Part of the strong observational evidence of the tidal strip¬ 
ping scenario includes: 


Carbon stars also in the MB area (e.g. Kunkel et al. 


2000 ). 


• Intermediate-age stellar populations in the Southern 
SMC ( Noel fc Gallart|2007 1. 

• Kinematically peculiar stars in the LMC that might 


be evidence of being captured from the SMC (Olsen et al. 


2011 ). 


• Possible evidence of extra-tidal stars at ~10.6° from the 
SMC centre ( Nidever et al.pOll ). 

• Old stars found in the MB ( Bagheri et al.|2013 ). 

• The red clump stars found in the OGLE IV field^ 
( Skowron et al.|2014 ). 

• The stellar structure in the eastern SMC identified by 


Nidever et al. (2013). These authors suggest that this is a 


stellar counterpart of the HI MB that was tidally stripped 
from the SMC around 200 Myr ago (during a close encounter 

^ http://ogle.astrouw.edu.pl/ 


with the LMC). This implies an increment in the self-lensing 
of the SMC stars, in line with the tidal origin for the mi- 
crolensing events reported towards the LMC by|Besla et al.| 

(Msl). 


The intermediate-age stars in the MAGIC I (Noel et al. 


2013a) in a held closer to the SMC make up 28% of all stars 


in that region. These are not present in helds in a direc¬ 
tion pointing away from the LMC. This provides potential 
evidence that these intermediate-age stars could have been 
tidally stripped from the SMC. 


With the advent of the new STElj^survey (Ripepi et al. 


2014) that will cover the main body of the SMC (32 deg^;, 
the MB (30 deg^) and a small part of the Magellanic Stream 
(2 deg^), evidence will likely continue to mount. 

In spite of the above striking evidence, there still re¬ 
main some results at odds with the tidal scenario. IDemersI 
|fc Battinelli| ( |1998| ) carried out the hrst extensive study of 
the Inter-Cloud stellar population Ending only young stars 


(~10-25 Myr old). A similar result was also found by Harris 


(2007) who studied twelve Helds in the MB Ending only old 
stars east of a~03h 18’ and (5~-74° and only young stars 
in their Inter-Cloud fields along the MB’s main ridge-line, 
with no evidence of an intermediate-age component. As they 


did not End evidence of tidally stripped stars, Harris (2007) 


concluded that the material stripped from the MCs into the 
MB must be almost purely gaseous, while any young stars 
must have formed in-situ in the MB. 

With the purpose of shedding light on the stellar con¬ 
tent of the MCs’ Inter-Cloud area, we have initiated the 
MAGIC (MAGellanic Inter-Cloud) project, aimed at disen¬ 
tangling the population age and distribution in the Inter- 
Cloud region as well as probing the faint outskirts of the 
MCs in the direction of the Bridge. As stated above, in 


MAGIC I (Noel et al. 2013a) we performed a quantita¬ 


tive Color-Magnitude Diagram (CMD) fitting analysis of 
two large Inter-Cloud fields, finding that some ~28% of 
the stars in one of the fields have ages between ~1 and 
10 Gyr old (intermediate-age population). [Harris ] ( [2007| ) 
likely missed this intermediate-age population because their 
CMDs were substantially shallower than the MAGIC CMDs 
that reach the oldest MS turnoffs, of key importance for 


spotting intermediate-age stars (see e.g. Noel 2008). The 
results from MAGIC I provide more evidence that these 
intermediate-age stars could have been stripped from the 
SMC, supporting the tidal scenario. 

In the present paper we focus on the stellar content of 
ten fields strategically located perpendicular to the MAGIC 
I field where intermediate-age population was found repre¬ 
sented by a purple triangle in Figure (field B2, located 
at ~8 Kpc from the SMC centre in the Inter-Cloud region). 
Five of the fields are radially away from the MB towards the 
North and five are radially away from the MB towards the 
South of field B2 in MAGIC 1. This arrangement allows us to 
differentially compare stellar populations inside the bridge 
region with directly neighbouring fields, allowing us to sep¬ 
arate distant SMC stars from those possibly stripped from 
the SMC’s interior. Our goal here is to probe the outskirts 
of the SMC in the MB direction. 


2 STEP: The SMC in Time: Evolution of a Prototype interacting 
late-type dwarf galaxy’ 
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This paper is organized as follows. In section we de¬ 
scribe the observations, data reduction and photometry. The 
analysis and results are discussed in section Finally, in 
section]^ we discuss our results and present our conclusions. 


2 OBSERVATIONS, DATA REDUCTION AND 

PHOTOMETRY 

Following our experience with previous observations for 
MAGIC I ( |Noel et al.|[M3a| , we obtained B and R band 
images of ten fields between the nights of 18*’’ to 26*’’ De¬ 
cember 2011 using the Wide Field Imager (WFI) on the 
2.2m telescope in La Silla. The WFI consists of 8 chips cov¬ 
ering a field of view of approximately 30’x30’ and has a 
pixel scale of 0.24”/pixel. Each B-band exposure was 1200 
seconds with 900 seconds for the R-band. An overview of the 
field positions can be seen in Figure with the MAGIC II 
fields depicted as squares showing the Northern and South¬ 
ern extensions from our base position. The Northern arm ex¬ 
tends radially away from the SMC, while the Southern arm 
extends perpendicular from the Magellanic Bridge (we will 
return to this Figure later). In Tablewe present the field 
name, RA/DEC centres of the observed fields, observation 
date, seeing estimate, filter, airmass, and the 50% complete¬ 
ness level of the data. N represent the fields located North 
of B2 and S those located South of B2, the increasing num¬ 
bers in S (STS5) represent fields of increasing Declination 
from the bridge ridge-line while the increasing numbers in N 
(N1-N5) represent fields of decreasing Declination. The over¬ 
laps between the MAGIC I field (B2) and the first Northern 
field is 16.8 arcminutes and 14.3 arcminutes for the Southern 
field. Each subsequent field has an overlap of approximately 
2-3 arcminutes. For the overlap with MAGIC I this corre¬ 
sponds to at 2300 stars and for the overlaps between fields 
we obtain an average of 780 stars to determine the photo¬ 
metric offsets. 

The data reduction steps performed here followed the 
same procedure as in MAGIC I. The basic reduction was 
performed by the Cambridge Astronomical Survey Unit 
pipeline (CASU; [Irwin fc Lewis|[2001| with nightly master 
bias and master flat frames from adjacent nights were com¬ 
bined to build rolling-average master frames for each night 
of data. In the past the ESO/WFI instrument has had prob¬ 
lems related to scattered light resulting in the need to apply 
a zeropoint correction map to the data, as per [Monelli et| 
jal.j ( |2003| ). With the installation of the baffle in 2006, this 
has dramatically reduced the scattered light incident on the 
camera. Attempts to measure the scattered light in WFI in 
2012, using a specially designed survey 086.A-9021(A), did 
not recover a scattered light signature as seen by previous 
authors. For this reason, we have not applied any correc¬ 
tion of this manner. The CASU pipeline was then used to 
generate aperture photometry and to determine the astro- 
metric solution by cross-matching with the 2MASS point 
source catalogue. As per MAGIC I, the astrometric solution 
is generally better than 0”.2. The final photometric cata¬ 
logue and artificial star tests were undertaken usi ng point- 
spread function photometry pipeline DOLPHOlj^ (jPolphiri 


^ http://purcell.as.arizona.edu/dolphot/ 


20001. DOLPHOT extracts the sources and determines the 


magnitudes for the two bands, B and R, simultaneously. This 
was run on all of the fields creating a catalogue of sources 
per field in both filters. The output files of dolphot include 
detailed information for each star including the position and 
global solution such as: Xi sharpness, roundness, and object 
type. For our purposes, we selected a value for y ^ 3.5, 
reasonable for good stars within relatively crowded fields, a 
value of (sharpness)'^^O.l for a perfectly-fit star, and object 
type (which ranges from 1 to 5) of 1 indicating a ‘good star’. 
These parameters allow us to remove spurious objects such 
as extended objects and cosmic rays from the photometry. 
Typically this corresponds to ~3000 objects removed per 
field from our area of interest in CMD space. 


Artificial star tests were used to determine the pho¬ 
tometric completeness of the observations and consisted of 
~48300 fake stars per chip, dolphot places each artificial 
star iteratively into the science image and then attempts to 
extract each them along with the real stars. Each success 
or failure of this process allows us to build up a profile of 
the photometric completeness. Following the procedure in 
MAGIC I, we fitted the ratio of the detected fake stars to 
total number of fake stars per magnitude bin with a Logistic 
function and from this we were able to determine the 50% 
completeness level of each frame. This information is then 
utilized in the CMD fitting routine MATCH which will be 
discussed in a later section. 


To calibrate the MAGIC II fields, we use the photo¬ 
metric solution for MAGIC I and from the overlaps between 
each successive field determine the correction for both the 
Northern and the Southern Arms. In this process, we tie 
the instrumental magnitudes of all fields in each arm to the 
N1 and SI field respectively, as this only requires a sim¬ 
ple offsets between fields. From this solution we then apply 
the remaining colour terms and offsets determined from the 
overlap with MAGIC I to correct the photometry to the 
same Johnson-Cousins system as MAGIC 1. We utilised the 
Markov Chain Monte Carlo code emcee ijForeman-Macke^ 


et al. 20131 to simultaneously calculate the best offsets in 


both B and R for all fields when determining their relation 
to N1 and SI. Similarly, when determining the colour term 
and offset between MAGIC I and Nl/Sl we also used em¬ 
cee. In both cases, we typically used over 1000 walkers with 
more than 5000 steps each. The solutions for the fits can 
be found in Table The errors in the offsets, as determined 
using emcee, are passed on in quadrature to each successive 
field. The evolution of the photometric error with magnitude 
for both the Northern and Southern Arms are shown in Fig¬ 
ure The spread in errors is greater for the B band as it has 
a much larger colour term than the R band, which is essen¬ 
tially zero. This has the effect of broadening the error profile. 
As expected though, at the bright end, the minimum error 
for each field is greater than the previous field. The resulting 
CMDs are shown in Figure [^and|^ where the CMDs corre¬ 
sponding to the Northern and Southern fields (with respect 
to B2 of MAGIC I) are presented. The red dashed lines in 
Figureshow the region used for our analysis here. In order 
to avoid the contaminants, such as background galaxies, at 
the faint end of the magnitude limit we only took the stars 
up to R=22.3. 
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Table 1. Summary of Observations. The zero-points of the data are 25.25 in B and 24.68 in R. The atmospheric 
extinction parameters, = 0.271 and = 0.07, were taken from < |Giraud et al.|2006| l and l |Tug|1977| |, respectively, 
as this gave the best B-R colour in comparison with established catalogues. The colour term applied to all fields 
was 0.223di o'ool ^ band and -O.OlSdi q'ooi 


Field RA DEC Date Seeing Filter Airmass 50% Mean E(B-V) Offsets 

Name Observed Completeness 


N1 

35.668180 

-71.99990 

2012-01-18 

1.01" 

E 

a] 



35.667640 

-72.00006 

2012-01-18 

0.97" 

F 

®[ 

N2 

36.166910 

-71.49983 

2012-01-19 

1.09" 

B 




36.166580 

-71.49979 

2012-01-19 

1.06" 

R 

N3 

36.500400 

-71.00002 

2012-01-20 

1.14" 

B 



36.500700 

-70.99997 

2012-01-20 

1.03" 

R 

N/ 

c 

d 

37.001040 

-70.49991 

2012-01-22 

1.42" 

B 




37.000260 

-70.50009 

2012-01-22 

1.20" 

R 

N5 

37.500900 

-70.00009 

2012-01-25 

0.91" 

B 



37.500630 

-70.00077 

2012-01-25 

1.05" 

R 

SI 


35.668100 

-73.50038 

2012-01-20 

1.07" 

B 



35.667680 

-73.50001 

2012-01-20 

1.02" 

R 

S2 


36.167360 

-73.99998 

2012-01-21 

CO 

00 

B 



36.168150 

-73.99988 

2012-01-21 

CO 

O 

R 

S3 


36.501480 

-74.49991 

2012-01-23 

1.65" 

B 



36.500620 

-74.50014 

2012-01-22 

1.16" 

R 

S4 


37.000630 

-74.99995 

2012-01-24 

1.14" 

B 



37.001710 

-75.00065 

2012-01-24 

1.06" 

R 

S5 


37.500590 

-75.49980 

2012-01-26 

1.05" 

B 



37.498980 

-75.50018 

2012-01-26 

O 

to 

R 


1.44 

25.06 

0.024 

-0.593 

H- 

OO 

bb 

OO 

-TO 

1.51 

24.28 


-0.724 

+ 0.012 

1.63 

24.78 

0.028 

0.296 

- 0.004 

1.76 

23.78 


0.119 

+ 0.004 

- 0.001 

1.59 

24.61 

0.033 

-0.038 

+ 0.003 

- 0.003 

1.73 

23.89 


0.062 

+ 0.003 

- 0.001 

1.40 

24.67 

0.036 

-0 01 

u-ui _o.o 03 

1.49 

23.95 


0.004 

+ 0.003 

- 0.004 

1.41 

24.89 

0.043 

0.101 

+ 0.003 

- 0.003 

1.49 

23.78 


0.082 

+ 0.004 

- 0.003 

1.55 

24.67 

0.040 

-0.605 

0.005 

0.021 

1.49 

23.89 


-0.778 

0.004 

0.01 

1.64 

24.39 

0.043 

-0.049 

+ 0.003 

- 0.003 

1.75 

23.23 


-0.068 

+ 0.002 

- 0.002 

1.71 

24.00 

0.043 

-0.051 

+ 0.002 

- 0.003 

1.63 

23.56 


-0.111 

+ 0.001 

- 0.001 

1.51 

24.84 

0.036 

-0.158 

+ 0.004 

- 0.001 

1.57 

23.89 


-0.021 

+ 0.002 

- 0.003 

1.52 

24.84 

0.037 

-0.104 

+ 0.002 

- 0.001 

1.58 

24.17 


0.005 

+ 0.001 

- 0.003 


BB#B/123_ES0878 
^ BB#Rc/162_ES0844 
^ N4, Chip 7, B offset: -0.5 l°:o03 
^ N4, Chip 7, R offset: 0.04 


2.1 CMD fitting 


The CMD fitting technique was carried out using the soft¬ 
ware package MATCH ( Dolphin|2002 1 as explained in Noel 
|et al.| ( |2013a| where the reader should refer for more details. 
In summary, the observed CMDs are converted into Hess di¬ 
agrams and compared with synthetic CMDs of model popu¬ 
lations from |Girardi et ah] ( |2004[ ). Theoretical isochrones are 
then convolved with a model of the completeness and pho¬ 
tometric accuracy in order to create the synthetic CMDs. 
We also obtained a foreground estimation based on a Galac¬ 
tic structure model given by MATCH that is included in 
the software as an extra model population. The software 
uses a maximum-likelihood technique to hnd the best linear 
combination of population models, resulting in an estimate 
of the star formation history (SFH) and the age-metallicity 
relation (AMR). In order to obtain the SFHs for the Inter- 
Cloud helds analyzed here, we used the colour range -0.5 ^ 
(B-R) 1.2 and the magnitude range R ^ 22.3, as marked 

by the red dashed lines in Figurej^ In this way, we avoid the 
inclusion of red low mass foreground stars in the fitting. To 
further test the influence of unidentified background galax¬ 
ies contaminating our sample in the region between 21 R 
^ 22.3, we calculated the ratio of objects in the red clump 
to the number of objects at the lower edge of the CMD (21 
^ R ^ 22.3). We find these ratios to be 10-12%. This is 
comparable with the ratio of ~11% for the wing SMC held 


from Noel et al. (20071 (qjOlll, see grey small diamond in 
Figure [^, where the contamination is known to be very 
low. This, together with the constraints given by dolphot 
described above, assures us that any background contami¬ 
nation has been minimised for our dataset. 

We recovered the SFH and metallicity for all the helds 
using age bins in the range from 10 Myr to 13 Gyr and metal- 
licities in the range -2.4:^ P^/H] 0 in bins of width 0.2 dex. 

We used a Salpeter initial mass function (IMF; since these 
helds contain predominately intermediate- and old-age stars, 
the choice of IMF is not critical as shown in |de Jong et al. 


20101 and assumed a 30% binary fraction (see 


Noel et al. 


20071. The htting software then hnds a linear combination 


of model CMDs that best ht the observed CMD. The result¬ 
ing relative star formation rates (SFRs) and the [Fe/H] as 
a function of time for the Southern and Northern helds are 
presented in Figures and [^respectively. 


3 RESULTS 

3.1 The star formation histories in the MAGIC II 
fields 

The upper panels of Figures and [^ show the SFHs for all 
helds. The SFR for each bin is the average of the values 
given by the hts at the different distances and the error bars 
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represent the complete range of SFR values found. As in the 
case of the MAGIC I fields, we find here a conspicuous old 
population (older than ~10 Gyr old) present in all fields. 
It is striking to note that the bulk of intermediate-age and 
young(er) stars increase as we move closer to field B2 of 
MAGIC I located in the vicinity of the MB main ridge-line. 
Very little star formation took place in the past 10 Gyr in 
the Southernmost and Northernmost fields. These results 
agree with the proposed numerical scenario where the MB 
is the result of a recent tidal interaction between the MCs 
and the MW ( Muller fc Bekki||2007 K The amount of HI (as 
seen from Figure [Ij) is similar in all the Northern arm fields 
as in the Southernmost field S5. 

The lower panels show the AMRs. The metallicities 
are only plotted for age bins with significant star forma¬ 
tion (SFR isO.2). The metallicity of each bin is the average 
of the fits weighted by the SFH. Error bars in the AMR 
represent the standard deviation. Horizontal bars in both 
the SFH and the AMR plots represent the width of the age 
bins used. The AMRs are consistent with being constant 
on average at [Fe/H]~-0.8 in agreement with the results of 
MAGIC I ([Noel et al.||2013a|) and with the Calcium triplet 


metallicities obtained for the outskirts of the SMC (Car- 
rera et al.|2008 1 and - within the quoted uncertainties - the 


LMC (Carrera et al. 20111. Recently, Bobbie et al. (20141 
found similar results using a factor 10 of more stars. Note 
that these AMRs obtained from CMD fitting should only be 
taken as indicative of the metallicity range. 

As shown in Noel et al. (2013aI, valuable quantitative 
information on the stellar distribution is given by the cumu¬ 
lative SFH that provides the fraction of stellar mass formed 
prior to a given time. In Figures and we show these 
cumulative SFHs for our 10 MAGIC H fields. The upper- 
j/-axis shows the redshift as a function of tim^ with the 
blue dashed lines mark the time when 50% of the stars were 
formed. Although this varies across the fields, most of the 
helds had 50% of their stellar population before 10-12 Gyr 
ago. 

The previously introduced Figure gives a visual pic¬ 
ture that shows the density HI map using data from the 
LAB survey of Galactic HI ( Kalberlaet al.|2005 l. The colour 
contours denote the HI emission levels integrated over the 
velocity range 80<v [km/s] <400. Each contour represents 
the HI column density twice as large as the neighbour con¬ 
tour (Jan Skowron, private communication). Our MAGIC H 
helds are overlapped and depicted as bright magenta, pink 
and off-white squares (these differences are explained be¬ 
low). MAGIC I held B2 is also shown in Figu re [l] a s a pur¬ 
ple triangle. |Harris| ( |2007[ ) helds, |Bica et al.| ( |2008[ ) associ¬ 
ations and clusters in the relevant right ascension (2h-4h), 
and OGLE IV helds 102 up to 125 from |Skowron et al.|2014] 
(without numbers for clarity purposes) are also overlapped. 


3.2 Southern Arm fields 

Here we concentrate on disentangling the age distribution 
in the Southern arm helds S1-S5 as seen in Figure The 
MAGIC H helds in bright magenta indicate that we hnd a 


The redshifts were obtained using the following cosmology: 
H(}=71, S1m= 0.270, f2A=0.730 (http://map.gsfc.nasa.gov/) 


fair amount of intermediate-age stars there (between ~27% 
and ~48%). Fields in pink show substantial amount of 
intermediate-age population (between ~18% and ~22%) 
while oh-white helds indicate that very little intermediate- 
age population is found there. The areas with larger amounts 
of intermediate-age stars coincide with the higher HI density 
column but they are more spread out than the HI. In these 
regions, there are a greater amount of young stars, stellar 
clusters and associations (Bica et al. 20081. As we move away 
from the higher HI areas, the amount of intermediate-age 
population diminishes until there are very few intermediate- 
age stars in the Southernmost held S5. The latter are some¬ 
how puzzling since there are young associations found in 
those parts of the MB. However, our results are in agreement 
with helds 110, 113 and 114 from Skowron et al. ( 2014[ ) (the 
OGLE IV helds overlapping our Southern Arm helds) who 
do not hnd young or intermediate-age component in those 
areas. Note that the MAGIC I held B2 presents some ~28% 
of intermediate-age stars. 

From the cumulative SFHs shown in Figure we hnd 
that helds SI and S2 show a signihcant amount of stars 
formed between 10 and 1 Gyr ago. The remaining stars 
formed in the last Gyr with a conspicuous increment (20% 
of the stars) in the past few Myr. This recent burst is the 
only clear star burst associated with a pericentric passage 
between the Clouds (~ 0.2 Gyr ago; Kallivayalil et al.|2006 |. 
The arrows in SI and S2 in Figureshow the last two peri- 
centric passage between the LMC and the SMC :^2.5 Gyr 
ago and ~250 Myr ago, respectively, from |Diaz fc Bekki| 
(2012). In their best model, Diaz & Bekki (2012) hnd that 
the SMC disk remains intact until the dynamical coupling 
of the MCs ~2.4 Gyr ago, at which point the Magellanic 
Stream and the Leading Arm are violently torn away by the 
strong tidal forces of the LMC. They claim that a second 
tidal encounter occurring ~250 Myr ago is the responsible 
for pulling the MB from the SMC disc. From our qualita¬ 
tive comparison, there seems to be a correlation between 
this last pericentric passage between the MCs but no clear 
connection can be made between the other enhancements in 
the SFR and the pericentric encounters between the MCs, in 
agreement with the previous results from Noel et al. (20091 
and MAGIC 1. 


3.3 Northern Arm fields 

We focus here in the age distribution in the Northern Arm 
helds N1-N5 from Figures and As in the case of the 
Southern Arm helds, there is a bulk of intermediate-age stars 
in these helds. From Figure |^it is clear that almost no stars 
younger than 1 Gyr old were formed, ruling out the in-situ 
formation for the Northern Arm. 

The density column of HI in the Northern Arm helds 
is considerable lower than in the Southern Arm. However 
we still see intermediate-age stars indicating that the young 
stars follow the HI gas but the intermediate-age population 
is considerably more spread out. 


3.4 Disentangling the Age Spread 

We now turn our attention to the correlation between young, 
intermediate and old-age stars in the MAGIC H helds and 
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the HI gas. Without velocity information, we cannot deter¬ 
mine conclusively if such correlations exist, but marginalis¬ 
ing over the unknown stellar velocities we can ask if such 
correlations are at least supported by the data. In Figure]^ 
we look for evidence for a correlation between the HI gas 
and young stars in our helds. We sum up the number of 
young MS stars in each pixel of the SMC gas map from 


Muller et al. (20031. Stepping through each velocity slice of 


these radio data, we multiplied the star and gas densities 
and summed all positive pixels. The sum of this convolution 
as a function of the velocity slice is shown in Figure We 
find three velocity peaks where the stars and gas density 
correlate strongly, suggesting that it is at least possible for 
the young stars and HI gas to be highly correlated. 

Unfortunately, we cannot repeat the above analysis for 
our intermediate age and old-stars since these cannot be so 
easily identified by their colour and magnitude alone. In¬ 
stead, we consider how the intermediate-age and old stars 
relate spatially to our young star population at the spatial 
resolution of each MAGIC-H held. This is shown in Figure 
|10[ The greyscale contour shows the percentage of young 
(left), intermediate (middle) and old-age (right) stars in each 
of our MAGIC II helds. These are overlaid on gas density 
contours taken from fSkowron et al. ( 2014[ ) as in Figure 

Notice that the young stars are clustered strongly 
around our MAGIC I held B2 (marked by the purple trian¬ 
gle in Figure]^. By contrast, the intermediate-age stars are 
more extended to the north and south, while the old stars 
are spread with almost constant density across all helds. 
Thus, if the young stars correlate strongly with the HI gas 
then the intermediate and old-age stars do not. We discuss 
the implications of this in section 

3.5 Surface Brightness Profiles 

We derived the surface brightness for the ten helds analyzed 
here using the selection box shown in Figure as red dashed 
lines. The surface brightness is measured in B-band as a 
function of distance from the mass centre of the SMC. We 
integrated the hux in the CMD subtracting the hux contri¬ 
bution from foreground MW halo giants and metal poor MW 
dwarfs and from background galaxies. In order to estimate 
and remove the contamination, we used the predictions from 
the TRILEG AL Galaxy model ( Girardi et al.|2005 l scaling 
the number of predicted MW stars to the number observed 
in our helds in those parts of the CMDs that were clearly 
devoid of MB Inter-Cloud stars (see |Gallart et al.|2004 and 
Noel fc Gallart|2007 for more details). 

We obtained two separate huxes: the total hux from 
stars older than 10 Gyr old (considered the old population) 
and the total hux from intermediate-age stars (between 1 
and 10 Gyr old). The results are shown in Figure [TT] where 
the squares represent the Southern helds and the circles de¬ 
pict the Northern helds. Open (red in the colour version) 
symbols denote old populations while hlled (green in the 
colour version) ones denote the intermediate-age popula¬ 
tions. Error bars represent the diherence between the values 
of the foreground hrrx in the hducial CMD area as predicted 
by the TRILEGAL code and after scaling by the ratio of 
observed and predicted huxes (see |Noel fc Gallart|2007 1. As 
seen from the hgure, both prohles (from the old and from 
the intermediate-age populations) are relatively hat within 


the errors instead of the expected exponential at inner radii 
from the SMG centre. This is in agreement with the density 
prohle from Nidever et al. ( 2011| who found that at these 
remote distances the star counts prohle deviates from the 
exponential ht, showing a “break” in the population that 
extends up to 10.6° from the SMC center. [Nidever et al.| 
(20111 argue that the origin of this break population is ei¬ 
ther a tidal debris or a “classical” bound stellar halo similar 
to the one around the MW and M31. We cannot, however, 
directly compare their analysis with ours since [Nidever et ah] 
(20111 deliberately removed their outer helds in the direc¬ 


tion toward the LMC. 

The most plausible explanation is that the close peri- 
centre passage between the LMC and the SMC ~200 Myr 


ago (e.g. Besla et al. 2007 Noel et al. 20091, believed to 


have stripped HI to form the MB, could have also tidally 
stripped stars of all ages from the SMC as indicated by the 
similar spatial distribution of the old and intermediate-age 
prohles in Figure The fact that intermediate-age stars 
show a brighter surface brightness prohle could mean that 
the number of these intermediate-age stars is larger than the 
old one or that the intermediate-age stars are brighter than 
the old stars. 

Finally, it is important to note that the intermediate- 
age and old stars in Figure pTj are similarly distributed, while 
in Figure[^they have clearly different spatial distributions. 
This apparent contradiction arises because in Figure we 
see that the young, intermediate and old age stars are dis¬ 
tributed diherently with respect to the HI gas in the Magel¬ 
lanic Bridge. By contrast, in Figure El we consider instead 
how they are distributed with respect to distance from the 
SMC. Over the small radial range of 6.5 to 8.5 Kpc from the 
SMC centre probed by our MAGIC H fields, we are unable 
to detect a difference in the radial profile of the intermediate 
and old-age populations (within our quoted uncertainties). 
It is likely, however, that over large radial ranges and/or 
with more complete spatial coverage we would see rather 
different distributions for the intermediate-age and old-age 
stars, as we do in Figure El 


4 DISCUSSION AND CONCLUSIONS 

We have presented the stellar content and distribution of a 
transversal strip perpendicular to the main HI ridge of the 
Tnter-Cloud’ region. The fields were chosen deliberately to 
allow us to trace the stellar populations of the Inter-Cloud 
radially out from the SMC since thus the contaminating halo 
Galactic stars can be removed and allowing us to probe the 
outskirts of the SMC in the MB direction. 

We confirm here the results of MAGIC I, finding that 
the stellar populations and mean metallicities of the fields 
are similar to fields to the South and West of the SMC but 
at smaller radii from the SMC center ( ~ 2.5 Kpc as in the 
case of Noel et al.|2009 1. The similarity between the MAGIC 
I and MAGIC II fields to some SMC fields that lie much 
closer in projection to the SMC is an indication that these 
MB fields contain tidally stripped SMC stars. 

In good agreement with MAGIC I and previous stud¬ 
ies, we found significant intermediate-age stars in the bridge 
region and young stars of a similar age to the last peri- 
centre passage in-between the clouds ~200-250 Myr ago, as 
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determined from their proper motions and radial velocities 
(see Figure [^. These recent enhancements are likely associ¬ 
ated with the recent pericentric passage between the clouds 
I Kallivayalil et ^|2006 l. However, it is interesting to note 
that we find no such prominent bump in the star formation 
at earlier times. This could owe to the difficulty of determin¬ 
ing orbits backwards in time | Lux et al.|2010 1, or the poorer 


temporal resolution of the SFH backwards in time (Dolphin 


2012), or could indicate that the two clouds have only just 
interacted for the very first time ( Diaz Bekki|2012 ). 

We showed here for the first time that the young, 
intermediate-age and old stars have distinct spatial distribu¬ 
tions with respect to the HI gas in the MB. Our key results 
are as follows: 


(i) The young stars correlate well with the HI gas suggest¬ 
ing that they were either recently stripped from the SMC 
or formed in-situ. Unfortunately, at our current CMD-age 
resolution these two options are degenerate. Based on their 
proper motions, the clouds likely interacted just ~200 Myr 
ago (Besla et al. 2012[ | which is comparable to the recent 
star formation burst (see Figure]^. Thus, these stars could 
have formed in the body of the SMC and had time to be 
tidally torn into the MB (along with the bridge gas); or 
they could have formed in-situ. We note that both possibili¬ 
ties could act in tandem, while both lead to young stars that 
are highly correlated with dense gas. 

(ii) The intermediate-age stars (1-10 Gyr old) are located 
mainly in the bridge region but are more spread out than 
the young stars, as seen in Figure 

(iii) The intermediate-age population has very similar 
properties to stars located ~2 Kpc from the SMC centre, 
suggesting they were tidally stripped from this region. 

(iv) Finally, the old stars extend to some ~9 Kpc from 
the SMC supporting the idea that all galaxies have a large 
extended metal poor stellar halo. 


The work presented here and in MAGIC I represent a 
complex picture for the MB region. There appears to be 
young stars that correlate with HI gas; a spatially distinct 
intermediate-age component; and a possibly smoother old- 
age component. It is important to emphasize that this study 
shows the power of deep resolved CMDs as a tool for em¬ 
pirically unravelling the past history of even rather complex 
galactic systems like the Magellanic Clouds. 

In a future paper, we will use stellar kinematics in or¬ 
der to pin down whether the gas and young stars perfectly 
correlate; to explore the role of ram-pressure on the gas; to 
test whether the young stars formed in-situ; and to further 
probe the likely tidally-stripped intermediate-age stars. 
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MAGIC I, B1 and B2, depicted as the two purple triangles. Fields from |Harris|2007l are also overlapped with the corresponding numbers 
as well as OGLE IV fields 106 up to 125 from [Skowron et al.|20f4| (we avoided the numbers here for clarity purposes). |Bica et al.|2008| 
clusters and associations in the area around our MAGIC II fields are represented with blue dots. The different shades in the colours of 
the MAGIC II fields describe the amount of intermediate-age stars present in each field: the off-white squares depict fields with little 
intermediate-age population, the light pink squares represent those areas where there is substantial intermediate-age stars, and the bright 
magenta squares where there is significant intermediate-age population. The colour contours represent the HI emission integrated over 
the velocity range 80< v [km/s] <400, where each contour represent the HI density column taken from the Leiden/Argentine/Bonn 
-LAB- survey of Galactic HI jKalberla et al.|2005[|; see Figure 8 from|Skowron et al.|201^for more details. 
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Figure 2. Magnitude versus Errors in Magnitude for all fields in MAGIC II. The left panels represent the Northern Arm and the right 
panels represent the Southern Arm. The top panels are from the B-band photometry and the lower panels from the R-band photometry. 
The vertical dashed lines delineate the approximate region used in the CMD fitting analysis. The fields (in the bright region) are ordered 
by increasing error from 1 to 5 in both the North and South. Alternatively the fields 1 through 5 are represented by the colours Blue, 
Green, Red, Cyan and Magenta. 
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Figure 3. CMDs of the Southern arm fields presented in order of increasing Declination from top left. Red dashed lines represent the 
region used for the SFH analysis. Isochrones were overlapped showing different evolutionary phases. 
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Figure 4. CMDs of the Northern arm fields presented in order of decreasing Declination from top left. Red dashed lines represent the 
region used for the SFH analysis. Isochrones were overlapped showing different evolutionary phases (see figure^for a visual clarification). 
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Figure 5. SFH and AMR for the five Southern MAGIC II fields. Clockwise from top left panel in order of increasing Declination SFHs 
and AMRs for the ten MAGIC II fields (see figure Each Figure shows the average SFRs over the total age range (upper panel) and 
the SFR-weighted metallicity, i.e. the AMR (lower panel). The vertical error bars in the upper panels give the full range of SFRs found 
in the fits while the vertical error bars in the lower panels give the standard deviation in the metallicities of the fits. Horizontal bars in 
both the upper and lower panels indicate the width of the age bins. 
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Figure 6. Same as figure[^for the fields North of B2 in MAGIC I (see figure^for a visual clarification). 
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Figure 7. Cumulative SFHs. The blue dashed horizontal line shows when 50% of the stars were formed. The upper y axis shows redshift 
for our current cosmology. The (tentative) pericentre passages that indicate possible past interactions between the LMC and SMC as 
taken from[^az Sz Bekki|201^are marked by arrows in SI and S2. 
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Figure 8. Same as figure[^for the fields in the Northern arm (see figure^for a visual clarification). 
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Figure 9. Spatial distribution of the MAGIC II young MS stars are binned into the same spatial grid as the [Muller et ar]|2003| gas 
maps. The x-axis represents the different velocity slices as per the jMuller et al.|2003| maps and the y-axis is the sum of the stellar density 
map multiplied by the gas density at that velocity. If the gas follows stars then the highest correlation between the maps occurs at ~180 
Km/s, while there are two other velocity peaks at ~193 and ~225 Km/s. 
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Figure 10. The percentage of young (left), intermediate-age (middle), and old (right) stars for our MAGIC II fields (grey contours). 
The coloured contours show the HI gas map from Figure^ 
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Figure 11. B-band surface brightness profile of the MAGIC II fields presented here. The intermediate-age (filled green symbols) and 
old stars (open red symbols) in each field were split in order to quantify the contribution of the different populations. Circles represent 
the Northern Arm fields while squares denote the Southern Arm fields. See figure [T] for a visual location of the fields. 
























